Phonon Scattering and Internal Friction in Dielectric and Metallic Films at Low 

Temperatures 
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We have measured the heat conduction between 0.05 K and 1.0 K of high purity silicon wafers 
carrying on their polished faces thin dielectric films of e-beam amorphous Si, molecular beam epi- 
taxial (MBE) Si, e-beam polycrystalline CaF2, and MBE CaF2, and polycrystalline thin metallic 
films of e-beam Al, sputtered alloy Al 5056, e-beam Ti, and e-beam Cu. Using a Monte Carlo sim- 
ulation to analyze the conduction measurements, we have determined the phonon mean free path 
within the films, and found all of them to be much shorter even than in typical bulk amorphous 
solids, with no exceptions. We have also measured the internal friction of these films below 10 K 
and found, however, their internal friction at low temperatures strikingly close to that of amorphous 
solids, both in magnitude and in their temperature independence, with the exception of the MBE 
Si and alloy Al 5056, whose internal friction is even much smaller than that of amorphous solids. 
The internal friction results indicate the phonon scattering in these thin films is the same as, or 
even much less stronger than, in other amorphous solids, according to the Tunneling Model. Thus, 
we conclude that the heat conduction measurements do not support the picture that the lattice 
vibrations of these films are glasslike, as had been surmised earlier for thin metallic films, on the 
basis of low temperature internal friction measurements alone [Phys. Rev. B 59, 11767 (1999)]. At 
the least, the films must contain additional scattering centers which lead to the very small phonon 
mean free path. Most remarkably, the MBE Si shows the same strong scattering of thermal phonons 
as do the other films, while having the negligible internal friction expected for a perfect film. The 
disorder causing the strong scattering of the thermal phonons in this film is completely unknown. 
The non-glasslike phonon scattering phenomena observed here in thin dielectric and metallic films 
deserve further investigations. 



I. INTRODUCTION 



Structure and perfection of thin films on substrates 
are still poorly understood. The purpose of the present 
investigation is to show that thermal phonons with wave- 
lengths on the order of 100 nm can be used as very sensi- 
tive probes of their disorder. It will be shown that strong 
phonon scattering occurs in a large number of films, in- 
cluding silicon films produced on silicon substrates by 
molecular beam epitaxy (MBE) . The nature of this dis- 
order is, however, not understood. 

It has recently been shown that the internal friction 
of crystalline metal films below 10 K resembles that of 
amorphous solids both in magnitude and temperature in- 
dependence, the so-called internal friction plateau 0. A 
possible explanation was that crystalline metal films have 
the same density of tunneling states as amorphous solids. 
In disordered crystals, such states are called glass- like ex- 
citations @. According to the Tunneling Model (TM) 
||, both the low temperature thermal conductivity be- 
low 1 K and the internal friction plateau below 10 K are 
determined by the same quantity, the tunneling strength 
C: 



C 



Pj 2 



(i) 



where P is the uniform spectral density of the tunneling 
states, 7 is their coupling energy to phonons, p is the mass 
density, and vt is the transverse sound velocity. C deter- 
mines the internal friction plateau through relaxational 
scattering of the elastic wave, and the phonon thermal 
conductivity through resonant scattering of the thermal 
phonons. This quantitative connection between inter- 
nal friction and thermal conductivity has been proven in 
many cases for bulk amorphous solids and also for dis- 
ordered crystals, and constitutes a major proof of the 
validity of the TM, as reviewed in Refs. Q and ||. 

We have recently described a technique by which we 
can measure thermal phonon scattering in thin films on 
substrates Using this technique, we have verified 
the quantitative connection between internal friction and 
thermal conductivity successfully for amorphous silica 
films |(| and for crystalline silicon layers which had been 
disordered by ion implantation to the point of amorphiza- 
tion @. We will use this technique here for a comparison 
with the internal friction on a variety of dielectric and 
metallic films. If the lattice vibration of the films are 
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TABLE I. Preparation parameters of the thin films used for the heat conduction measure- 
ments in this work. The Cu film has a thin (100 A) adhesive layer of Ti between it and the 
silicon substrate. The alloy Al 5056 film contains, by weight, 5.2% Mg, 0.1% Mn, and 0.1% 
Cr. The films were deposited onto either one or both of the wide polished substrate faces 
(see Ref. B]), as indicated in the last column ("one" or "both"). 
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indeed glass-like, the phonon scattering should be de- 
termined by the same tunneling strength as determined 
from the internal friction. 

In addition to the phonon scattering in the metal 
films studied previously in internal friction we will 
also present measurements of both internal friction and 
phonon scattering in amorphous Si (a-Si) films, in dielec- 
tric crystalline CaF2 films, and in a crystalline Si film 
produced by MBE, which is expected to contain fewer 
defects than any of the other films. 

II. EXPERIMENTAL MATTERS 
A. Thin Films 

The thin films for heat conduction measurements were 
deposited either on Czochralski-grown, (111) oriented sil- 
icon substrate surfaces, or on float-zone refined, (100) ori- 
ented ones. For internal friction measurements, the thin 
films were deposited onto the double-paddle oscillators, 
which were float-zone refined, (100) oriented. All sub- 
strates were of high purity, and were double-side polished. 
Film thickness was determined by calibrated vibrations 
of a 6 MHz plano-convex quartz crystal. When possible, 
it was double-checked with a step surface profiler. De- 
tails on the films for heat conduction measurements are 
contained in Table I. 

In order to eliminate surface contaminations to ther- 
mal conduction measurements, samples were, when ap- 
propriate, either put through an RCA clean || or cleaned 
in a hot sulfuric acid solution [pi. For the internal fric- 
tion measurements, all films were deposited directly onto 



double-paddle oscillators after cleaning of the substrate 
by diluted HF solution, except for the MBE Si film. Be- 
cause of the stringent cleaning requirements prior to the 
MBE deposition in the UHV chamber, which were not 
suitable for double-paddle oscillators, the MBE Si film 
was deposited onto a wafer from which an oscillator was 
subsequently fabricated. The fabrication process involves 
heating the wafer to 850°C for 20 minutes. Thus, the 
MBE Si film for the internal friction measurement was 
considered as annealed. All the other annealing processes 
were done in the MOS area of the Cornell Nanofabri- 
cation Facility and were preceded by a stringent RCA 
cleaning, as described in Ref. pp| , in order to avoid any 
contamination of the silicon which is known to occur dur- 
ing annealing under regular clean laboratory conditions. 
This annealing process will be referred to in the following 
as "MOS-cleaned-and-annealed." 

Since we are primarily interested in the thermal 
phonon mean free path in the film from the heat con- 
duction measurement, we would like to minimize phonon 
scattering at the film-substrate interface and at the free 
surface of the film. As the film-substrate interface is ex- 
pected to be much smoother than the free surface of the 
film, we first consider the free surface roughness of the 
films studied in this work. Table II presents the root- 
mean-square (RMS) roughness of the free surfaces of the 
films studied as determined by atomic force microscopy 
(AFM). Table II also presents the dominant thermal 
phonon wavelength at 1 K based on the Debye speed 
of sound of the materials listed. The wave length of the 
thermal phonons is the relevant length. As can be seen, 
the RMS roughness is small compared to the length scale 
of the thermal phonons; and so, the free surfaces of the 
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TABLE II. Comparison of the free surface RMS 
roughness as determined by AFM measurements 
with the dominant thermal phonon wavelength Adom 
at 1 K based on the Debye speed of sound, vd, for 
the thin film samples studied. 
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a Taken from Ref. jllj or estimated using Eq. |^ with v t 
listed in Table III. 
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FIG. 1. Measured thermal conductivity of a high purity Si 
substrate (large faces polished and thin faces sandblasted, see 
Appendix) multiplied by its thickness, 300 iim: solid trian- 
gles; and the calculated electronic thermal conductivity of a 
Cu film (see text) multiplied by its thickness, 0.2 iim: solid 
line. At 50 mK, only about 5% of the heat is carried in the 
film (by electrons); at higher temperatures, the percentage 
drops as T~ 2 . 



films should have little effect on the overall scattering of 
thermal phonons and thus on the overall conclusions of 
this work. In section III C, we will show experimental ev- 
idence that phonon scattering at the surfaces/interfaces 
is indeed negligible relative to that in the films, by study- 
ing films thickness dependence. 

B. Methods 

In contrast to the conventional thermal conductivity 
measurements on bulk metals, our thermal method, as 
applied to metal films, has the advantage of solely de- 
termining the phonon heat transport in the form of a 
phonon mean free path in the film. In our investigation 
the metal films act primarily as phonon scatterers rather 
than as heat conductors because of their relatively small 
thickness to that of the substrate. At low temperatures, 
most of the heat in normal metals is carried by electrons. 
This thermal conductivity can be calculated using the 
Wiedemann- Franz-Lorenz law (reviewed in Ref. ||l2f ) and 
an appropriate electrical resistivity. Because heat trans- 
port is parallel to the film-substrate interface, the amount 
of heat carried by the film or the substrate depends on 
their relative thermal resistance, which is the inverse of 
thermal conductivity multiplied by length and divided by 
cross-sectional area. Since the film and substrate have 
the same length and width and differ only in thickness, a 
comparison of the products of thermal conductivity and 
thickness is enough to determine which carries most of 
the heat, which is shown in Fig. ^ using a 0.2 thick 
Cu film as an example. The low temperature thermal 
conductivity of the Cu film was determined with the 



Wiedemann-Franz-Lorenz law using a room temperature 
electrical resistivity of 1.6 x 10~ 6 Slcm and a residual 
resistivity ratio of 2, measured in our laboratory on sim- 
ilar films |l],|l^|. Fig. [l] demonstrates that the substrate 
phonons are the dominant heat carriers, primarily be- 
cause the thickness of the substrate is so much greater 
than that of the film. The same is true for the other 
metal films studied in this work, even more so when they 
become superconducting in the temperature range inves- 
tigated here. 

Thermal phonon mean free paths were determined, 
using a Monte Carlo (MC) simulation, from heat con- 
duction measurements between 0.05 and 1.0 K by the 
technique described in Ref. [||, denoted as ^fii m (HC) m 
the following. We mention briefly that additional scat- 
tering mechanisms such as scattering from free surface 
roughness can be included in the simulations, should 
that become necessary. For those who wish to determine 
^fiim(HC) horn heat conduction measurements below 1 K 
without having to resort to performing their own simula- 
tions, we provide information in the Appendix, using the 
results of our MC simulations on a film-substrate sample 
with dimensions as typically used in our work. 

We can also predict the phonon mean free path from 
internal friction measurements, denoted as 4im(TM)! if 
we assume that the film has the low energy excitations 
that are common in amorphous solids (and no other scat- 
tering centers) within the TM model. In the present 
work, the low-temperature internal friction of thin films 
is measured with double-paddle oscillators vibrating in 
their antisymmetric mode at ~ 5.5 kHz, which have ex- 
ceptionally small background damping as described pre- 
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viously 0|L5|]. Thin films increase the internal friction 
of the paddle oscillator, Q~l ddic - From this, the internal 
friction of the film, lm , is determined by |l4|] 



Q 



film 



Gsub^sub 
3Gfil m £film 



(Q- 1 



paddle 



(2) 



where t and G are thicknesses and shear moduli of sub- 
strate and film, respectively, and Q~ b is the internal fric- 
tion of the bare paddle (including the mounting losses). 
Gfii m is assumed to be equal to that of the bulk mate- 
rial The specific model used to obtain 4im(TM) 
from the internal friction of a film is the TM, originally 
proposed by Anderson, et al. [jiHf , and independently 
by Phillips fl7| , expanded for elastic measurements by 
Jackie |[8| . The TM connects the thermal phonon mean 
free path, £, with the internal friction plateau, Qq , as 
follows. From Ref. 
conductivity A is 



||, the expression for the thermal 



A : 



-C v v D £ 



2.66fc| 



(3) 



where, in the gas-kinetic picture, C v is the low tempera- 
ture specific heat per unit volume, vd is the Debye speed 
of sound, ks is Boltzmann's constant, h is Planck's con- 
stant, and T is the temperature. Note that 



Qo 1 — 



(4) 



where G is defined in Eq. [T]. Substituting for C v within 
the Debye model of the phonon spectrum 19 1, Eq. || be- 



comes 



Vi 



(1.59 x HT i2 [sK])— 



assuming the empirical relation flll| , [20|| : 

v t ~ 0.9u_d, 



(5) 



(G) 



where [sK] are units of second and Kelvin. These equa- 
tions provide the means to predict (within the TM) what 
^fiim(TM) should be if the internal friction plateau of the 
film, Qq 1 film , is known. To repeat, Eq. ^| assumes that 
the internal friction plateau is due to the presence of 
glassy states in the film, and that no defects other than 
the glass-like excitations scatter the thermal phonons. 
The validity of these assumptions will be tested for 
the films investigated here by comparing £fii m (HC) with 

^film(TM)- 



10- 



0.5/_tm 
a-Si 



10" 




O^DOOCOOOCOt 



0.5/im 

a— Si, annealed 



background 



0.4^m 



10- 1 10° 10 1 10 2 

temperature (K) 

FIG. 2. The internal friction of a bare double-paddle oscil- 
lator (solid curve "background" ) and of such oscillators carry- 
ing e-beam a-Si and MBE Si films. Note the negligible effect 
of the MBE film. The annealing of the e-beam a-Si film was 
done at 700° C for 1 hr under the MOS-cleaned-and-annealed 
condition (see text). The MBE Si film had been annealed at 
850° C for 20 min (see text). 
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FIG. 3. Internal friction of an e-beam a-Si film, before and 
after annealing, compared to that of bulk a-Si02 (solid curve). 
The bulk a-Si02 data, measured at 4.5 kHz, are taken from 
J.E. Van Cleve, Ph.D. thesis, Cornell, published in Ref. Q. 
The double-headed vertical arrow indicates the range of the 
temperature-independent internal friction plateau, measured 
on a wide range of bulk amorphous solids as reviewed in 
Ref. il. 



III. RESULTS AND DISCUSSION 

A. Silicon Films 

Since we are searching for tunneling states in thin films, 
we start with e-beam a-Si, a highly disordered film known 



to have such states |14j , and will compare it with a crys- 
talline silicon film produced by MBE which is expected 
to be a simple extension of the silicon lattice. Fig. 
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shows the internal friction of a bare double paddle oscil- 
lator, called "background," and of the same kind of oscil- 
lator carrying the films. As expected, the MBE film has 
negligible internal friction, while the e-beam a-Si films, 
both as-deposited and annealed, lead to a considerable 
increase of the internal friction. The MBE Si was mea- 
sured only in the annealed state as explained in Section 
II A. 

From the change of the internal friction of the paddle 
carrying the films, the internal friction of the films, Qguj,, 
can be determined using Eq. ^, and is compared to that of 
bulk a-SiC>2 in Fig. ||. The annealing of the e-beam film 
causes almost an order of magnitude reduction in Qgi m , 
while the Qgi m of the annealed MBE film is too small 
to be determined. Annealing at 700°C for 1 hr leads to 
almost complete crystallization of a 0.5 /jm thick e-beam 
a-Si film pl| . The internal friction confirmed that ~ 90% 
of the low energy excitations had been removed. Below 
1.0 K, however, the internal friction increased, indicative 
of a contamination in c-Si Jl0| , although the most strin- 
gent MOS-cleaned-and-annealed process, as described in 
Section II A, was strictly followed. Most probably, an 
impurity was trapped on the substrate surface during 
mounting in the e-beam evaporator which is located out- 
side the MOS area. This impurity, trapped underneath 
the a-Si film, survived the MOS cleaning and led to the 
contamination during the annealing. In contrast, there 
is no such problem for the MBE Si film, and hence no 
contamination is observed. 

Table III summarizes Qq 1 film , v t , 4im(TM)( T ) given by 
Eq. ^[ for the films presented here and below. This en- 

TABLE III. The internal friction plateau Qq film , 
the transverse speed of sound Vt, and the thermal 



phonon mean free path I 



film(TM) 



Eq. |^, where T is measured in Kelvin 



(T) predicted by 
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a Taken from Ref. 
b Taken from Ref. 
c Taken from Ref. 
d Taken from Ref. 
e There is no Qq 1 
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value for the MBE Si film because 
the internal friction of this film was not detectable (see 
Fig. ||), and hence ishnCTU) is expected to be very long. 



ables us to compare the internal friction and the phonon 
mean free path directly within the TM. 

Fig. U shows ^fiim(HC) m the MBE Si and e-beam a-Si 
films both in as-deposited and annealed states obtained 
from the heat conduction measurements. The dotted line 
represents ^fiim(TM) of the e-beam a-Si based on the in- 
ternal friction measurement (see Table III). For e-beam 
a-Si, £fii m (HC) is significantly smaller than 4i m (TM) which 
assumes that the scattering occurs by tunneling states 
alone. Evidently, £gWHC) cannot be used to test for the 
existence of glassy excitations in e-beam a-Si film. In 
addition to the tunneling states, other scattering centers 
must be present. The situation may be similar to e-beam 
a-SiC>2 films in which the additional thermal phonon scat- 
tering was explained by cracks or voids Q. It is well 
known that a-Si films produced by e-beam evaporation 
contain similar defects [^3|-^5|. After annealing, only a 
small increase of ^fii m (HC) above 0.3 K and even a de- 
crease below that temperature can be seen in Fig. [|. This 
change of ^awnc) upon annealing, which cannot be ex- 
plained within the TM, is interpreted as resulting from 
a combination of a decreased scattering by the tenfold 
smaller number of tunneling states, an increased scat- 
tering by the contaminants, and possibly a change in 
scattering from the defects of unknown nature which had 
been noticed already in the film prior to annealing. Obvi- 
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FIG. 4. Phonon mean free path 4im(HC) of as-deposited 
and annealed silicon films. As-deposited e-beam a-Si: solid 
circles; annealed (700° C, 1 hr) e-beam a-Si: Open circles; 
as-deposited MBE Si: solid squares; annealed MBE Si (500° C, 
1 hr followed by 700°C, 1.5 firs): Open squares. The dotted 
line is 4im(TM)(r), the TM prediction based on the internal 
friction plateau of the as-deposited e-beam a-Si. Dashed lines 
are guides for the eye. 
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ously, these results cannot be used to extract any knowl- 
edge about the annealing of these unknown defects. They 
do, however, provide further evidence for the presence of 
the contaminants introduced into crystalline silicon dur- 
ing annealing under all but the most stringent conditions, 
thus further emphasizing the need for their identification 
and control JlO[ . 

Very surprisingly, Fig. ^ shows that the phonon scat- 
tering in MBE Si film both before and after anneal- 
ing is very large as well. Since the contamination that 
plagued the e-beam a-Si films is not an issue here, sec- 
ondary ion mass spectroscopy (SIMS) was performed by 
R. Reedy at the National Renewable Energy Laboratory 
(NREL) in order to identify other possible chemical im- 
purities. The following chemical elements were detected: 
boron, < 10 16 cm -3 ; nitrogen, ~3x 10 16 cm" 3 ; carbon, 
~3x 10 17 cnr 3 ; and hydrogen, ~2x 10 18 cm" 3 . These 
concentrations were similar in the MBE film and the sili- 
con substrate. Only the oxygen contents differed between 
substrate (~5x 10 17 cm" 3 ) and film (~5x 10 18 cm" 3 ). 
For all these detected impurities, an anneal (500° C, 1 hr 
followed by 700° C, 1.5 hr) caused no measurable change 
of their concentrations. Since no evidence for such scat- 
tering was observed on bare silicon samples from the 
same batch as used in these experiments jjj], the only 
possible impurity scatterer is the oxygen. But the ten- 
fold increase of oxygen in the MBE film should not lead 
to a noticeable decrease of the experimental mean free 
path t (defined in the Appendix, Eq. [?]), and thus to a 
decrease of ^fiim(HC) given the relatively small thickness 
of the film, unless the oxygen in the film somehow acts as 
a much stronger scattering center. Phonon scattering in 
oxygen-doped silicon has been found to depend on heat 
treatment |^6|,^7j at frequencies in excess of ~ 300 GHz. 
But, phonons in this frequency range carry heat predom- 
inantly above 3 K, and no evidence exists for phonon 
scattering by oxygen at lower frequencies (corresponds 
to T < 1 K). 

In an attempt to detect any evidence for structural dis- 
order in this MBE film, an x-ray diffraction (XRD) anal- 
ysis was performed by M. Sardela and D. Cahill at the 
University of Illinois (Champaign-Urbana). High resolu- 
tion open-detector scans around the Si (004) peak showed 
no difference between measurements conducted on the 
film side and on the substrate side of the MBE Si film- 
substrate sample. Full width at half maximum values 
were found to be almost identical on both sides, and no 
diffuse scattering or any disorder feature on the film side 
was seen. In addition, a triple axis reciprocal space map 
around the Si(004) peak on the film side also could not 
detect any diffuse distribution nor asymmetry of the Si 
peak. These observations speak against the existence of 
grains with different orientations which might cause the 
phonon scattering. 

The near T _1 temperature dependence may be sug- 
gestive of scattering by sessile dislocations as reported, 



for example, by Wasserbach in plastically deformed bulk 
copper, niobium, and tantalum p^j. If we assume that 
the coupling between dislocations and phonons is simi- 
lar, the dislocation density in the MBE film would have 
to range between 10 10 and 10 13 cm -2 , which seems rather 
high for MBE silicon. At this point, no search for dislo- 
cations in this MBE film has been undertaken. 

MOS-clcaning-and-annealing of the MBE film as de- 
scribed above leads to an increase of ^mm(HC) shown in 
Fig. [|, although it remains below the mean free path pre- 
dicted even for e-beam a-Si, except at the lowest temper- 
ature. Thus, noticeable disorder, other than the tunnel- 
ing states, remains in the MBE Si film even after anneal- 
ing, although the internal friction of the annealed MBE 
Si film shown in Fig. ^ give no evidence for any low en- 
ergy excitations. Thus, the only firm conclusion we can 
draw at this point is that the defects in the MBE film 
are not glass-like. 

B. CaF 2 Films 

As was just shown, crystalline films produced by crys- 
tallizing an a-Si film or by MBE deposition show little or 
no evidence for tunneling states in low temperature inter- 
nal friction. It is therefore surprising that a 0.6 \m\ thick 
film of crystalline e-beam CaF2 increases the damping 
of the double paddle oscillator by more than one order 
of magnitude, see Fig. ||. The internal friction of the 
film itself, as compared with that of a-Si02 shown in 
Fig. |^, is nearly temperature independent and close to 
the range found for all amorphous solids studied to date 
(with the exception of certain hydrogenated a-Si films, as 
discussed in Ref. |ljl). Thus, the large internal friction 
observed previously in polycrystalline metal films 0] ap- 
parently also occurs in some crystalline dielectric films. 
Assuming that its cause is glass-like excitations, we can 
again predict an ^fii m (Twn for the CaF2 film, shown as 
the dotted line in Fig. [t|. The measured ^fii, n (HC) f° r 
an identical CaF2 film, also shown in Fig. (7), is more 
than two orders of magnitude smaller than predicted by 
the TM. We also measured the phonon mean free path 
of another crystalline CaF2 film which was prepared by 
MBE technique at 750° C substrate temperature to im- 
prove its structure. However, the internal friction of the 
MBE CaF2 film cannot be measured because of the dif- 
ficulty of preparing such a film on silicon paddle oscil- 
lators and meeting the requirements of special cleaning 
and fabrication at the same time. The ^fii m (HC) °f the 
MBE CaF2 film, though larger than that of the e-beam 
one, is still smaller than that predicted by the TM for 
the e-beam CaF2 film. Furthermore, it is smaller than 
that of an a-Si02 film prepared by wet-thermal oxida- 
tion, in which the structure is much improved in com- 
parison with e-beam a-SiC>2, and the phonon scatter- 
ing is determined solely by the glassy excitations. The 
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FIG. 5. The internal friction of a bare high purity silicon 
double-paddle oscillator (solid curve "background") and of 
such a paddle carrying the e-beam CaF2 film on the polished 
silicon surface. 




CaF 2 

aaaaaaaaaaaa a 



10 -5 ■■ 

10° 10 1 10 z 

temperature (K) 
FIG. 6. Internal friction of the e-beam CaF2 film compared 
to that of bulk a-Si02 (solid curve). The bulk a-Si02 data, 
measured at 4.5 kHz, is taken from J.E. Van Cleve, Ph.D. the- 
sis, Cornell, published in Ref. Q. The double-headed verti- 
cal arrow indicates the range of the temperature-independent 
internal friction plateau, measured on a wide range of bulk 
amorphous solids as reviewed in Ref. |5|. 
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FIG. 7. Phonon mean free path ^fiim(HC) of two different 
CaF2 films. MBE CaF2i solid circles; e-beam CaF2: solid 
squares. The dotted lines are the TM prediction based on in- 
ternal friction measurements for e-beam CaF2 and for thermal 
a-Si02, respectively. Data for the thermal (open squares) and 
e-beam (solid stars) a-Si02 films are taken from Ref. M . For 
thermal a-Si02, good agreement is shown between -£m m (HC) 
and ^fli m (TM) , which had also been found in ion-implanted sil- 
icon [Ql, as mentioned in Section I. Dashed lines are guides 
for the eye. 



though probably more highly ordered than an e-beam 
one, contains disorder because the pseudomorphic epi- 
taxial growth is known to break down at film thickness 
exceeding 10 nm [^9|,^o|, leading to structural relaxation. 
Internal stresses are also expected to result from differ- 
ential thermal contraction as the sample is cooled from 
the deposition temperature. There is, however, no con- 
vincing evidence for glass-like excitations in this CaF2 
film. As in the e-beam a-Si film, some unknown scatter- 
ing process masks the effect of the glass-like excitations, 
if they exist at all. 



C. Metal Films 



^fiim(HC) °f the thermal a-SiC>2 film agrees perfectly with 
that of the TM's prediction, just as one would expect 
for bulk a-SiC>2, as shown in Fig. ^ (see Ref. ^ for 
details). The ^mm(HC) of the MBE CaF 2 film locates 
between those of the thermal a-SiC>2 film and a macro- 
scopically as well as microscopically disordered e-beam 
a-SiC>2 film, in which the thermal phonon scattering 
is not dominated by the tunneling states, also shown 
in Fig. 0. We suggest that this MBE CaF 2 film, al- 



Fig. H shows ^ n i m (HC) f° r three e-beam Al films, 0.2, 
0.4, and 0.6 /jm thick. The absence of any significant 
dependence on the film thickness validates the assump- 
tion used in our analysis that the scattering occurs pre- 
dominantly within the films and not at the interfaces, 
an assumption which so far had been based only on the 
smoothness observed on the free surfaces as listed in Ta- 
ble II. The dotted line is the prediction for 4im(TM) based 
on the internal friction of the e-beam Al film reported 
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temperature (K) 
FIG. 8. Phonon mean free path 4im(HC) f° r e-beam Al 
films that are 0.2 itm (open triangles), 0.4 rim (open squares), 
and 0.6 /im (open diamonds) thick. The dotted line is the TM 
prediction based on internal friction (see Table III). Dashed 
lines are guides for the eye. 
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FIG. 9. The thermal conductivity of 5% deformed bulk 
Al. The thermal conductivity of a bulk a-Si02, taken from 
Ref. p2| . along with the phonon thermal conductivity of the 
0.2 thick e-beam Al film converted from their phonon 
mean free path shown in Fig. @, is shown for comparison. 



earlier [JjJ (see also Table III). As for the two previous ex- 
amples, the observed phonon scattering far exceeds the 
scattering expected on the basis of the TM. 

In Ref. pi], it had been shown that the low tempera- 
ture internal friction of an Al film on a Si substrate was 
very similar to that of heavily deformed bulk Al. It was 
therefore interesting to compare the phonon mean free 
path in the film with that observed in deformed bulk Al. 
For that purpose, a 99.999% pure polycrystalline Al rod 
(2.5 mm in diameter and 25.7 mm long) was first an- 
nealed at 560°C and subsequently stretched by 5%. Its 
thermal conductivity, measured by the standard steady- 
state technique, is shown in Fig. || along with that of bulk 
a-Si02- The steep rise of the thermal conductivity of the 
bulk Al above 0.1 K is caused by the onset of heat trans- 
port by normal state electrons. However, below that tem- 
perature, heat is expected to be carried predominantly 
by the lattice, and a temperature dependence similar to 
that of the bulk glass (a-Si02) is observed. Although the 
magnitude is three times smaller, it still falls within the 
glassy range in thermal conductivity, see Fig. 1 in ref. 
plf . The phonon thermal conductivity of the 0.2 /mi e- 
beam Al film, as calculated from ^fii m (HC) m Fig- Hi is 
also shown in Fig. l| Above 0.1 K, the thermal conduc- 
tivity of the deformed bulk Al and the phonon thermal 
conductivity of the e-beam Al film show the difference be- 
tween the heat transport by electrons and phonons, and 
by the phonons alone, separated here experimentally for 
the first time. Below 0.1 K, the phonon thermal conduc- 



tivity of the e-beam Al film is very close to that of the de- 
formed bulk sample. This suggests that the defects which 
scatter the phonons in the film are very similar to those in 
the heavily deformed bulk sample. The same conclusion 
had been reached previously in internal friction measure- 
ments as stated above. The defects causing the internal 
friction had been tentatively identified as dislocations or 
dislocation kinks Q. It is tempting to suggest that the 
thermal phonons in the films are scattered by the same 
defects. Since we see in Fig. || that ^fii m (HC) an d ^fiim(TM) 
of the e-beam Al films are not connected by the TM, we 
can conclude that the same holds for the deformed bulk 
Al because of the similarity in the internal friction and 
phonon mean free path between the thin films and the 
bulk samples. Thus, the non-glasslike phonon scatter- 
ing phenomena observed in this work are not limited to 
thin films alone. In addition, the defects or the mecha- 
nisms causing the resonant scattering of thermal phonons 
in heat conduction and those leading to the relaxational 
process in internal friction may not even be related, as 
shown by the following observation. 

The alloy Al 5056 in bulk form has an excep- 
tionally small low temperature internal friction, even 
as a sputtered film (Table III), which has been ex- 
plained by dislocation pinning The ^fii m (HC) 
in this film, however, is still close to that of all 
other metal films, see Fig. It follows that pin- 
ning of dislocations has no influence on the thermal 



8 



1000 





a. 



s 
I 



100 



10 




0.1 1 
temperature (K) 

FIG. 10. Phonon mean free path ^fu m (HC) f° r a 0.4 /im thick 
e-beam Al film: open squares; for a 0.5 /jm thick sputtered 
alloy Al 5056 film: solid circles; for a 0.1 ^m thick e-beam 
Ti film: solid triangles; and for a 0.1 /im thick e-beam Cu 
film: solid squares. The thermal conductivity of for the 5% 
deformed bulk Al is converted to its mean free path: open cir- 
cles. The labelled dotted lines are TM predictions of 4i m (TM) 
based on internal friction measurements (see Table III). The 
solid curve is the phonon mean free path in bulk a-Si02 taken 
from Ref. |B2|. Dashed lines are guides for the eye. 



phonon scattering. We conclude that the mechanisms 
causing the internal friction and the thermal phonon scat- 
tering are not understood. 

A comparison between ^fii m (HC) an d 4im(TM) is shown 
in Fig. [lO] for four metal films: Al, alloy Al 5056, Ti, 
and Cu, with the phonon mean free path of a-Si02 for 
comparison. The apparent lack of correlation between 
^fiim(HC) an d 4im(TM) enable us to generalize the same 
conclusion from the Al films to other metallic films, which 
is that if glass-like lattice vibrations exist in them, their 
effect is masked by the unknown defects. 



As observed in internal friction Q, 



film(HC) 



is unaf- 



fected by superconductivity (T c is 0.4 K for Ti, 0.92 for 
alloy Al 5056 Q, and 1.2 K for Al, ). It is concluded that 
phonon scattering by electrons is unimportant. Klemens 
has derived an expression for the phonon-electron scat- 
tering coefficient P in terms of the electron-phonon scat- 
tering coefficient E |34|] . Using the value for E measured 
by Berman and MacDonald for pure copper |35| , |36| ], we 
calculate £a\m (of phonons being scattered by electrons) 
at 1 K to be ~ 15 fim. This phonon scattering rate (due 
to electrons) is more than an order of magnitude less than 
the phonon scattering rate observed in Fig. ^ for the Cu 
film. Note that the calculation of 15 /im should not be 



taken too seriously as its assumptions |3j of the adia- 
batic principle, of a phonon Debye spectrum, and of a 
free electron gas may not be adequate at these tempera- 
tures for a thin polycrystalline Cu film with a residual re- 
sistivity ratio of 2 [fij. Nevertheless, this estimate agrees 
with our observation that electron-phonon interaction is 
not significant in our experiment. 



IV. CONCLUSIONS 



Measurements of the thermal phonon mean free path 
on films of amorphous and MBE Si, of polycrystalline 
and MBE CaF 2 , of pure metallic Al, Cu, and Ti, and 
of the metallic alloy Al 5056 below 1.0 K have revealed, 
in all cases, similar strong phonon scattering. Scatter- 
ing by surface and interface roughness can be excluded, 
since nearly the same ^ mm (HC) has been observed in Al 
films of different thicknesses. In searching for the origin 
of this phonon scattering, we have also measured the low 
temperature internal friction of the Si and CaF 2 films 
(that of the metal films had been measured previously, 
Ref. Q) and also the thermal conductivity of a bulk Al 
rod after a 5% plastic elongation. In all cases, £fii m (HC) 
was found to be much smaller than 4im(TM) based on 
the internal friction and assuming that the lattice vi- 
brations are glass-like. The discrepancy is particularly 
striking for the MBE Si film in which no internal fric- 
tion was observed, yet £a\m(RC) was similar to that found 
in all other films. In this case, the phonon scattering 
is particularly puzzling since the film is expected to be 
structurally more perfect. In all other films, macroscopic 
defects like grain boundaries, voids, cracks, or disloca- 
tions may be the cause for the phonon scattering. In the 
deformed bulk Al, the phonon mean free path was found 
to be equal to that in thin Al films. Since in the bulk 
sample, individual dislocations or aggregates thereof are 
likely phonon scatterers, they may also be the cause for 
the scattering in the films. However, dislocation motion, 
presumably tunneling, which has been invoked to explain 
the internal friction of deformed Al and of Al films (see 
Ref. UH) is an unlikely cause for the thermal phonon 
scattering since the same ^fii m (HC) m Al was also found 
in the alloy Al 5056, in which dislocation motion appears 
to be suppressed, resulting in a greatly reduced internal 
friction. In conclusion, both internal friction and phonon 
scattering have been shown to be sensitive probes for thin 
film disorder, including that in MBE Si. The nature of 
such disorder and the mechanisms by which it affects the 
elastic and thermal properties are completely unknown. 
No evidence for the existence of glass-like lattice vibra- 
tions has been detected. 
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Appendix 

The technique used in this investigation for the mea- 
surement of the thermal phonon mean free path in thin 
films has been described before ||. Although the exper- 
imental schematic resembles that of a thermal conduc- 
tivity measurement, see Fig. [ll], it should be emphasized 
that our experiment leads directly to a thermal phonon 
mean free path, rather than to a thermal conductivity 
A, from which the mean free path £ has to be calculated 
using the gas-kinetic expression 

A = ^C v vl, (7) 

which requires knowledge of the specific heat C v of the 
heat carrying excitations or phonons traveling with an 
average velocity v. In amorphous solids, for example, 
this C v cannot be measured. It can only be calculated 
from v through the use of the Debye model. 

The analysis of the heat conduction measurements 
on the silicon substrate carrying the film requires a 
Monte Carlo simulation which, though straightforward, 
is nonetheless time-consuming. By strictly adhering to 
the specifics as listed below (including sample and clamp 
geometry), one can extract the phonon mean free path 
in a film, £m m (called ^fii m (HC) m this paper), from £ exp 
(the experimentally measured phonon mean free path of 
the film-substrate sample) without having to repeat any 
MC simulations, as will be shown in this Appendix. 

Fig. [T^ is a plot of the results of MC simulations on a 
film-substrate sample with dimensions typically used in 
this investigation. For any particular simulation, £fum/x 
is an input parameter where x is the film thickness. The 
output parameter is ^mCi a simulated phonon mean free 
path of the film-substrate sample. To determine £mm, 
one sets an £ exp equal to an £mc in Fig- to find the 



silicon sampie 




FIG. 11. Schematic of 1-D heat conduction experiment; the 
sample is cleaved from a high purity commercial silicon wafer 
(orientation (111) or (100)) with both large faces polished; 
the thin faces are sandblasted as described in Ref. M. Also 
shown is a ballistic path of a thermal phonon from the silicon 
substrate through a thin film as modeled in the MC simula- 
tions. 
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FIG. 12. Typical plot of £mm/x versus £mc where x is the 
film thickness. Normalizing £m m with respect to x allows the 
same plot to be used for different film thicknesses as long as 
the sample and clamp geometry has not changed (see text for 
details). The actual code of the Monte Carlo programs may 
be found in Ref. §§. 

corresponding £mm/x] multiplying by x then yields the 
desired value, £f l \ mi that directly corresponds to the mea- 
sured 4x D - 
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In order to use Fig. [T^ as above, the film-substrate 
sample must be mounted as shown in Fig. O, Further- 
more, the four thin sides of the high purity silicon sub- 
strate must be completely roughened (by sandblasting, 
for example) while the two wide faces must be smooth 
from Syton polishing, for example), as explained in Ref. 

Any film of thickness x must cover both wide faces 
entirely and uniformly. The height of the sample above 
the top of the base clamp should be 44.5 mm, the width 
of the sample 7 mm, and the thickness of the substrate 
0.279 mm. The height of the heater, cold thermome- 
ter, and hot thermometer clamp above the top of the 
base clamp should be 41 mm, 4 mm, and 25 mm, respec- 
tively. The two interfaces between the heater clamp and 
the sample (areas of contact) should each be of dimen- 
sions 0.279 x 2 mm 2 while the four interfaces between 
the thermometer clamps and the sample should each be 
0.279 x 1.5 mm 2 . The effect of varying these details are 
discussed in Ref. ||. 
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